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Abstract 

Hydrodynamic fluctuations at non-zero temperature can cause slow relaxation toward equilib- 
rium even in observables which are not locally conserved. A classic example is the stress-stress 
correlator in a normal fluid, which, at zero wavenumber, behaves at large times as t _3//2 . A novel 
feature of the effective theory of hydrodynamic fluctuations in supersymmetric theories is the pres- 
ence of Grassmann-valued classical fields describing macroscopic supercharge density fluctuations. 
We show that hydrodynamic fluctuations in supersymmetric theories generate essentially the same 
long-time power-law tails in real-time correlation functions that are known in simple fluids. In 
particular, a i~ 3 / 2 long-time tail must exist in the stress-stress correlator of A^ = 4 supersymmetric 
Yang-Mills theory at non-zero temperature, regardless of the value of the coupling. Consequently, 
this feature of finite-temperature dynamics can provide an interesting test of the AdS/CFT cor- 
respondence. However, the coefficient of this long-time tail is suppressed by a factor of 1/Y?. 
On the gravitational side, this implies that these long-time tails are not present in the classical 
supergravity limit; they must instead be produced by one-loop gravitational fluctuations. 
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I. INTRODUCTION 



Holographic AdS/CFT duality implies that properties of strongly coupled large- N c quan- 
tum field theories can be deduced by doing calculations in classical (super) gravity [1-4]. The 
correspondence is believed to hold at non-zero temperatures, and can be helpful in extract- 
ing information about both supersymmetric and non-supersymmetric theories [5]. The most 
well-known example of AdS/CFT duality is A/"=4 supersymmetric Yang-Mills theory with 
gauge group SU(N C ), in 3 + 1 dimensions, which is believed to be dual to type II B string 
theory on the AdS$ x S* 5 background [4]. However, to date, there are very few physical 
properties which can be independently calculated on both sides of the duality and thus used 
as non-trivial tests of the finite temperature version of the correspondence. 

In this paper, we focus on the low-frequency real-time dynamics (rather than static ther- 
modynamics) of a finite temperature field theory. The relevant degrees of freedom are 
hydrodynamic fluctuations, by which one means those degrees of freedom whose relaxation 
time diverges with wavelength. 1 A suitable version of hydrodynamics is the appropriate form 
for an effective theory characterizing these degrees of freedom. As we will discuss in detail, 
one consequence of hydrodynamic fluctuations at non-zero temperature is the presence of 
long-time power-law tails in real-time correlation functions of conserved currents [6-8]. For 
correlations in the spatial parts of these currents at zero wavenumber, one finds t~ d l 2 behav- 
ior in d spatial dimensions, or equivalently, non-analytic cu d ^ 2 terms in the small-frequency 
behavior of the associated finite-temperature spectral densities. 2 These long-time tails are a 
generic feature of systems which behave as fluids on arbitrarily long time and large distance 
scales; their existence, and the value of the power-law exponent, are insensitive to the mi- 
croscopic details of the theory. Hence, long-time tails in real-time thermal correlators may 
be added to the small set of observables about which one can make firm predictions even in 
strongly coupled theories. 

We will apply our general results to the long-time behavior of thermal correlation func- 
tions in SU(N C ) J\f=4 supersymmetric Yang-Mills theory. Because this is a scale invariant 
quantum field theory, it cannot have phase transitions at any non-zero temperature (in infi- 
nite volume). 3 In particular, it must act like a fluid at all non-zero temperatures, and hence 
should have a valid hydrodynamic description of its long-time dynamics. Consequently, one 
should (for reasons we will discuss) expect the zero- wavenumber stress-stress correlator to 



1 In normal fluids, the relevant hydrodynamic degrees of freedom are fluctuations in energy and momentum 
density, or equivalently temperature and local fluid velocity. Hydrodynamic variables can also include 
fluctuations in charge densities associated with any globally conserved currents, as well as order parameter 
phase fluctuations (Goldstonc modes) in theories with spontaneously broken continuous symmetries. 

2 More precisely, the spectral density divided by ui, which is related to the power spectrum of thermal 
fluctuations, has a non-analytic \lo\^ 2 ^ 2 term at low frequency, unless d = 4ra + 2 in which case an 
additional ln|o;| is present. In many respects, long-time hydrodynamic tails are analogous to the non- 
analytic chiral logarithms which appear in one-loop chiral perturbation theory (at zero temperature) due 
to quantum fluctuations of massless Goldstone bosons. 

3 The 7V C — > oo theory has a phase transition in finite volume. This is irrelevant for our considerations. See 
Ref. [5] for more detailed discussion of the phase diagram. 
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show power-law relaxation, 

J d 3 x ({T^(t,x),T kl (0)}) ~r 3 / 2 , (1) 

as t — > oo. 

On the gravitational side of the AdS/CFT correspondence, the stress-stress correlator 
can be extracted from the absorption cross-section cr(u) for scattering of gravitons by non- 
extremal three-branes [9, 10]. This calculation was performed by Policastro and Starinets 
[11]. In the case when the graviton frequency uo was small, Policastro, Son and Starinets 
[12] interpreted the zero-frequency limit as the shear viscosity of M = 4 supersymmetric 
Yang-Mills plasma. In Ref. [11], the small frequency behavior of the cross-section was found 
to have the form a(uo) = <r(0) + 0(uo 2 ). However, if the stress-stress correlator does exhibit 
the power-law tail (1), then AdS/CFT duality implies that the corresponding graviton cross- 
section should be non-analytic at small u and behave as a(uo) — er(0) ~ M 1 / 2 - 

This apparent contradiction is the motivation for this paper. Our goal is to determine the 
coefficient of the long-time tail in the stress-stress correlator (1), and in analogous correlators 
involving spatial parts of other conserved currents. To do so, we will have to construct the 
correct effective description of low-frequency, long-distance dynamics in supersymmetric 
theories. We will find that supersymmetry (or superconformal) invariance of a quantum 
field theory does not preclude the existence of long-time tails, in accord with hydrodynamic 
expectations. However, in large N c gauge theories, we will find that a l/iV c 2 suppression 
appears in the amplitude of long-time tails. Thus, on the gravity side of the correspondence, 
the small-frequency non-analyticity can not be seen in classical supergravity (considered in 
Ref. [11]), but should emerge from a one- loop calculation on a plane-symmetric AdS black 
hole background. 

The paper is organized as follows. Section II reviews what is meant by an effective hydro- 
dynamic theory and summarizes the properties of hydrodynamic fluctuations in typical high 
temperature relativistic theories. Included is a discussion of why hydrodynamics predicts 
that fluctuations in the stress tensor T % \ at zero wavenumber, decay as t~ 3 / 2 at long times 
due to their coupling to both sound waves and transverse momentum density fluctuations. 
In Section III we construct effective hydrodynamics for supersymmetric theories, and argue 
that long-time tails are necessarily still present. The treatment of hydrodynamic fluctuations 
in sections II and III is applicable to general field theories. In Section IV we specialize the 
general results to M=A supersymmetric Yang-Mills theory in the context of the AdS / CFT 
correspondence. 

II. EFFECTIVE HYDRODYNAMICS 

Hydrodynamics may be viewed as an effective theory describing the dynamics of a ther- 
mal system on length and time scales which are large compared to any relevant microscopic 
scale [13, 14]. It is a classical field theory, whose fields can be regarded as expectation 
values in some non-equilibrium thermal ensemble of microscopic quantum operators aver- 
aged over spatial volumes large compared to microscopic length scales. The beauty of a 
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hydrodynamic description is that it applies to any system which acts like a fluid on suffi- 
ciently long distance scales, regardless of the strength of microscopic interactions. We will 
be concerned specifically with near-equilibrium behavior and the resulting dynamics of small 
perturbations about some equilibrium state. 



A. Hydrodynamic variables 



To construct the appropriate effective hydrodynamics for a particular physical system, 
one has to identify the relevant low-frequency degrees of freedom (hydrodynamic variables) 
and find the equations of motion which govern their dynamics. 4 As in any effective field 
theory, the symmetries, and symmetry realizations, of the underlying microscopic theory 
will constrain the form of the effective theory. All remaining dependence on the microscopic 
theory will be isolated in the particular values of some finite set of adjustable parameters 
appearing in the effective theory. In the case of hydrodynamics, these input parameters 
include transport coefficients (viscosity, diffusivity, conductivity, etc.) and equilibrium ther- 
modynamic functions. 

The degrees of freedom in a hydrodynamic description will include a minimal set of fields, 
whose thermal expectation values distinguish the possible equilibrium states of the theory. 
Among these are the energy and momentum densities, which we will denote as 

e = T 00 , vr^r , (2) 

as well as the charge densities n a = J° of any other conserved currents J^. 5 The argument 
for this is as follows. 

The usual grand canonical ensemble, with statistical density operator 6 

/3 = Z -l e /3(«,P-+^JV ) (3) 



built from the conserved charges 

P" = J d 3 x T M {x) , N a = J d 3 x J° a (x) , (4) 

describes a manifold of time-independent equilibrium states in which variations in the ther- 
modynamic parameters j3, u u , and fi a produce space-independent variations in the energy 
density, momentum density, and other charge densities. In other words, infinite wavelength 



4 Because the effective hydrodynamic theory describes dissipative dynamics, it is easier to work directly 
with the equations of motion than with a classical Lagrangian formulation. 

5 Until otherwise stated, we assume that such symmetry currents are ordinary bosonic vector fields (i.e., 
not supercurrents). 

6 We use a ( — h++) metric convention. (3 = 1/T is the inverse temperature, u v is the rest frame 4-velocity 
(satisfying u 2 = —1), and \i a are chemical potentials. We assume that translation invariance is a symmetry 
of the theory, and that this symmetry is not spontaneously broken. We work in three (flat, infinite) spatial 
dimensions throughout our analysis, but all results generalize trivially to d > 3 spatial dimensions. 
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variations in these densities have infinite relaxation time (precisely because they are den- 
sities of conserved charges). But in any local, causal theory, this implies that arbitrarily 
long wavelength variations in these densities must have relaxation times which diverge with 
wavelength. Consequently, fluctuations in conserved charge densities remain relevant vari- 
ables on arbitrarily long time scales, and must be retained in an effective hydrodynamic 
theory. 7 

If the theory under consideration is in a phase with spontaneously broken continuous 
symmetry, then the phase (or orientation) of the relevant order parameter is also needed to 
uniquely characterize equilibrium states, and fluctuations in the order parameter orientation 
will become another hydrodynamic degree of freedom. 8 Well known examples include su- 
perfluid helium, where the phase of the condensate wave function becomes a hydrodynamic 
variable and is responsible for the appearance of second sound [14]. Another example is the 
chiral limit of QCD, whose hydrodynamic variables include the SU(Nf) orientation of the 
chiral condensate, fluctuations in which describe pions [16]. At a second order phase transi- 
tion, long wavelength fluctuations in the magnitude of an order parameter acquire divergent 
relaxation times (due to critical slowing down), and also become relevant hydrodynamic de- 
grees of freedom. Finally, for gauge theories in a Coulomb phase, long wavelength magnetic 
fields have divergent relaxation times and must be retained in a hydrodynamic description; 
the effective theory in this case is termed magneto-hydrodynamics. 

In the following discussion we will assume, for simplicity, that all these complications are 
absent. That is, we assume that the theory under consideration is in a phase of unbroken 
global symmetry, is not in a Coulomb phase, and is not sitting precisely at a second order 
phase transition. Also for simplicity, we will assume that the theory possesses charge conju- 
gation symmetries under which any global conserved charges N a transform non-trivially. We 
further assume that the equal-time (or more generally, imaginary-time) thermal correlation 
functions of the theory exhibit a finite correlation length. This is basically just a restatement 
of our assumed unbroken global symmetry, absence of Coulomb phase gauge fields, and non- 
critical behavior. Finally, we assume that the theory under consideration is an interacting 
theory which describes a sensible equilibrating thermodynamic system — no discussion of 
hydrodynamic behavior is applicable to a non-interacting theory! All these assumptions 
hold in typical field theories without U(l) gauge fields, at sufficiently high temperatures. 9 



7 For theories with approximate symmetries, an approximately conserved charge density will only act like 
a hydrodynamic degree of freedom on time scales which are short compared to the mean time between 
charge non-conserving reactions. 

8 Despite the pervasive misuse of the phrase "spontaneously broken gauge symmetry" , Higgs phases of 
gauge theories are not examples of spontaneous breaking of any physical symmetry [15], and do not 
imply the existence of an enlarged manifold of physical equilibrium states, or the presence of additional 
hydrodynamic degrees of freedom. 

9 This includes non-Abelian gauge theories with most any matter field content, because interacting scalar 
fields, fermions, and non-Abelian gauge fields all develop finite spatial correlation lengths at high temper- 
ature. 
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B. Constitutive relations 



The resulting hydro dynamic description for this general class of theories takes the form 
of exact local conservation laws for the conserved currents, 

d,T^(x) = 0, (5) 
drJ£{x) = 0, (6) 

together with constitutive relations, valid on sufficiently long time and distance scales, ex- 
pressing the fluxes of conserved quantities (i.e., spatial parts of conserved currents) as local 
functionals of the hydrodynamic variables themselves (the densities of conserved charges). 10 
As in any effective theory, possible terms appearing in these constitutive relations must 
be consistent with the symmetries of the underlying theory, and may be classified accord- 
ing to the number of spatial derivatives, and powers of fields representing departures from 
equilibrium. This is the analogue of the usual power-counting in an effective field theory. 

We will focus on the dynamics of fluctuations away from some charge conjugation invari- 
ant equilibrium state, for which all chemical potentials vanish. Therefore, the charge density 
n a (x) will only be non- vanishing due to some perturbation away from the equilibrium state 
of interest. Similarly, if we work in the rest frame of the system, then the momentum density 
tt(x) will also be non- vanishing only due to a departure from our given equilibrium state. 

The only terms, linear in deviations from our equilibrium state, which can appear in the 
resulting constitutive equation for the charge fluxes (i.e., the spatial parts of the conserved 
currents Jg) have the form 

j a = -D ab Vn b , (7) 

where D = ||-D a &|| is, in general, a matrix of diffusion coefficients characterizing, in a linear 
response approximation, the flux induced by a spatially non-uniform charge density. These 
diffusion coefficients are input parameters to the effective hydrodynamics. 11 Terms with 
three or more spatial derivatives are also allowed by symmetry but, for sufficiently long- 
wavelength fluctuations, will be negligible compared to the above single derivative terms. It 
should be emphasized that the constitutive relation (7) is not an operator identity, unlike 
Eq. (6). Rather it is a relation between expectation values of operators in non-equilibrium 



More generally, one may formulate constitutive relations which, in addition to terms involving hydrody- 
namic variables, also contain noise terms representing the influence of short-distance degrees of freedom 
on the quantity of interest. Such noise terms convert the hydrodynamic equations of motion into Langevin 
equations, and allow the effective hydrodynamic theory to generate the correct equal time correlations of 
long wavelength equilibrium fluctuations. However, these noise terms will not be relevant for our purposes, 
and will be ignored throughout our discussion. 

In terms of equilibrium correlation functions in the underlying microscopic theory, diffusion coeffi- 
cients are given by the Kubo formula (Dx)ab — \ lim^^o / d 4 x e lu}x ° {\{Ja(x), J{(0)}), where (•••) = 
Z~ 1 tr[e~^ H ■ ■ ■] is an equilibrium thermal average, \ = ||Xa&|| is the charge susceptibility matrix, 
Xab = dn a /d((3iib) = [(N a Nb) — (N a )(Nb)] /V, and V is the spatial volume. The Kubo formula makes 
manifest the Onsager relation stating that D\ is a symmetric matrix. 
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states, valid only when those operators are averaged over a volume large compared to relevant 
microscopic scales. 

Of course, terms can also appear in the constitutive equation for j a which are non-linear 
in the departure from the given equilibrium state. At quadratic order 12 it is possible to add 
a term of the right symmetry which involves no spatial gradients, and is proportional to the 
product of charge and momentum density fluctuations, 

ja = ~D ab Vn b + k n a 7T . (8) 

Since this new term involves no spatial derivatives, the coefficient k is completely determined 
by thermodynamic derivatives which stay within the manifold of equilibrium states. Namely, 
an equilibrium state with non-vanishing chemical potentials, when viewed in an arbitrary 
reference frame, will have an energy-momentum tensor of the perfect fluid form, 

T i»> = ( e + uV + Vg»\ (9) 

and conserved currents 

J» = n a u^ (10) 

where n a , e, and V are the equilibrium charge densities, energy density, and pressure, re- 
spectively, in the rest-frame of the fluid, and is the rest-frame 4- velocity. 13 Hence, for 
infinitesimal boosts away from the rest frame, the momentum density and charge flux are 
related to the flow velocity via 

7r = (e + V)v, j a = n a v, (11) 

implying that the coefficient k appearing in the constitutive relation (8) equals (e+V)' 1 . For 
later convenience, we will denote the enthalpy density, which is the sum of energy density 
and pressure, as 

w = e + V, (12) 

so k = w~ l and 7r = w v. 

In the constitutive relation for the spatial stress T tj , terms without spatial derivatives 
are also completely determined by equilibrium thermodynamics. Let /3, e, and V denote the 
inverse temperature, energy density, and pressure, respectively, of the chosen equilibrium 
state with vanishing charge and momentum densities. Then the stress tensor of a nearby 
equilibrium state with energy density e + 5e, momentum density 7r, and charge densities n a , 
expanded to second order in deviations away from the reference state, has the form 

V + v 2 s 5e + U (fe) 2 + \ Z ab n a n b - v, 



lequil. 



.2 77 



S W 



+ — , 13 

w 



Terms higher than quadratic are irrelevant for our purposes. This will be discussed further at the end of 
section HE. 

The local rest frame of any thermal system, at a particular event x, may always be defined as the frame 
in which the momentum density 7r(x) vanishes. The local flow velocity in any other frame is then defined 
as the velocity needed to boost to the local rest frame. This implies that the local 4- velocity u v (x) is, in 
general, the timclike eigenvector of T^ v {x). 
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where Se = T 00 — e is the fluctuation in energy density, v s = (dp ' /de) 1 / 2 is the speed of 
sound, £ = d 2 P/de 2 , and S ab = d 2 P / dn a dn b . (Energy derivatives are to be evaluated at 
constant charge density, and vice- versa.) No terms linear in n a can be present due to the 
assumed charge conjugation invariance of the reference equilibrium state. An exercise in 
thermodynamic derivatives shows that 

2 dP dV fdsY 1 

v s=m=W\W = 1 C ' (14) 

i = [(l + 2v 2 )C- 1 -wT HHH ]/P, (15) 
We have defined C and x as the mean square fluctuation in energy or charge per unit volume, 



C = (# 2 ) con n /V, (17) 
Xab = (N a N b ) conn /V, (18) 

and T HHH and Y NNH as "amputated" third order connected correlators of energy and charge, 



= <# 3 > conn /V, (19) 

^ b NH = C" 1 (x-'U {X- I )w (N a , N b ,H) co jV , (20) 



with V the spatial volume, (...) denoting an expectation in the reference equilibrium state, 
and x _1 the matrix inverse of x = llXabll- Note that the heat capacity per unit volume 
c v = ds/dT = /3 2 C, while the charge susceptibility dn a /dfi b = (3 Xab- 

Terms linear in momentum density can also appear in the constitutive relation for stress, 
but only when combined with a spatial gradient so as to yield a rank-2 tensor. Two linearly 
independent structures are possible, proportional to the shear or divergence of the vector 
field, so that 



IT 2 

V + v 2 s 5e + \i (Se) 2 + \ E ab n a n b - v 2 — 



+ 



7T 7T 



W 



- 7c S ij V • 7T - lr] (V* TT j + V j 7T l - | 5 ij V • 7r) . (21) 

The coefficients 7^ and 7^ are conventionally written as 

7c = i, 11 = 1, (22) 
w w 

with ( and i] the bulk and shear viscosities, respectively. 14 

All other terms, consistent with symmetries, which could be added to the constitutive 
relations (8) or (21) involve either more gradients, or more powers of fluctuations away from 



Kubo formulas for the bulk and shear viscosities are £ = lr lim^^o J d*x e tbJt (^{V(x), V(0)}), and rj — 
I, lim^o / d 4 xe wx ° (±{s ij {x), s ij (0)}), where V = \T il and s ij = T^-V6 ij arc the trace and traceless 
parts of the stress tensor, respectively. 
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equilibrium relative to the terms included in (8) or (21), and so are negligible for sufficiently 
small, long- wavelength fluctuations. 15 

The above constitutive relation for the stress tensor simplifies considerably in scale- 
invariant theories. The tracelessness of the energy- momentum tensor, = 0, is an operator 
identity in such theories, valid in both equilibrium and non-equilibrium states (at vanishing 
chemical potentials). Consequently, in any scale invariant theory the equation of state is 
exactly e = 3V, the speed of sound v s = and the thermodynamic curvatures £ and 

S vanish identically, 16 as does the bulk viscosity ( . Therefore, the constitutive relation for 
the stress tensor (21) in a scale invariant theory takes the simpler form 

2*3 = gij (p + I Se) + 4 (vr* n j - 1 5 ij 7T 2 ) - 7r , (V 1 n j + V 3 7r* - 1 5 l] V • tt) . (23) 



C. Linearized Hydrodynamics 



If one retains only terms linear in fluctuations in the constitutive relations (8) and (21), 
then it is trivial to solve the resulting linearized hydrodynamic equations. For charge density 
fluctuations, inserting j = —D'Vn into the conservation law d t n — — V ■ j yields the 
diffusion equation (dt — DV 2 )n = 0. (We have suppressed indices distinguishing multiple 
charge densities, but in general D is a matrix of diffusion constants.) Inserting a spatial 
Fourier transform, 

n(t,x) = j ^ e *-n(f,fe), (24) 

immediately gives 

n(t,k) = e~ k2Dt n(0,k), (25) 

showing that fluctuations with wavevector k relax diffusively at a rate k 2 D which vanishes 
quadratically as k — > 0. Momentum density fluctuations may be decomposed into longitu- 
dinal and transverse parts, 

/d 3 k r - i 

^-^e lfe - a! [fe7r ll (t,fe) + 7r x (t,fc)J, (26) 

where k ■ n ± (t,k) = 0. For transverse momentum fluctuations, the conservation relation 
d t n % = — Vj T y plus the constitutive relation (21) again combine to give a diffusion equation, 
(dt — 7»;V 2 ) 7r ± = 0, so that 

n ± (t,k) =e- fc27 "*7r ± (0,fc). (27) 



It turns out that the non-linear interactions of hydrodynamic fluctuations cause the coefficients of higher 
order terms to become scale dependent quantities, just like couplings in a typical effective field theory 
[17]. 

Hence, for scale invariant theories w = § e, C = 4e/(3, T HHH = jg/3/e 2 , and T NNH = x^/w. These 
relations also follow by simple thermodynamics, given a free energy F ~ VT 4 and susceptibility Xab ~ T 3 . 



9 



Because energy and momentum densities are both conserved quantities, fluctuations in en- 
ergy and longitudinal momentum density are coupled. One finds that both quantities satisfy 
the damped oscillator equation (— d 2 — 7 S k 2 dt — v 2 k 2 ) Se(t, k) = 0, where 7 S = 7^ + §7^, 
and that \k\ n n (t, k) = id t 5e(t, k). In the long wavelength limit (7 s |fc| <C v s ), one has weakly 
damped sound waves, 

\v s 5e(t,k) J 6 V -sm(|fc|u 8 t) cos(\k\v a t) J \v s 6e(0,k) J ' 1 ' 

propagating at the sound speed v s , and whose energy attenuates at a rate of k 2r y s . 



D. Real-time correlators 



The non-equilibrium linear response results (25), (27) and (28), describing the relaxation 
of specific perturbations away from equilibrium, may be converted into equivalent results for 
real-time correlation functions characterizing the spectrum of fluctuations in the equilibrium 
thermal ensemble. 17 Consider, for example, the two-point correlation of charge density 
fluctuations, 18 

A na n b (t, k) = J d 3 x e~ ik x (± {n a (t, x),n b (0, 0)}) . (29) 
For t > 0, one may insert the linear response result (25) and find 19 

A nn (t,k) = \\A nanb (t,k)\\=e- k2Dt X (k), (30) 
where x(fc) = ||x a &(fc)|| is the variance matrix of equal-time charge density fluctuations, 

Xab (k) = J d 3 xe~ ik x (i{n a (x),n 6 (0)}>. (31) 

In the small-A; regime where a hydrodynamic description is valid, one may replace this vari- 
ance by its small k limit, which is the charge susceptibility matrix x introduced previously, 

lim xab(k) = [(N a N b ) - (N a ) (N b )} /V = Xah ■ (32) 
fc— >o 



This is sometimes referred to as Onsager's postulate. 

The effective hydrodynamic theory, in which fields are classical, cannot distinguish between the sym- 
metrized correlator (|{n(t), n(0)}) and the Wightman correlator (n(t)n(O)). In frequency space, this 
difference is smaller than the correlators themselves by a factor of Huj/T, which is negligible in the low 
frequency domain where the hydrodynamic description is valid. We have chosen to write expressions 
involving the symmetrized correlator just because this correlation function is always real (for Hermitian 
operators). 

For t < 0, one may use time-reversal (or CPT plus rotation invariance) to show that A nn {— t, k) = 
A nn (t, fc). Consequently, a Fourier transform in time gives A nn (uj 7 fc) = 2k 2 D [uo 2 + (fc 2 !?) 2 ]^ 1 \ ■ Com- 
bined with the Ward identity d M {\{ J»(x), J"(0)}) = 0, one finds associated hydrodynamic forms for 
the charge-current correlator, A nJ i (ui, fc) = 2ui k l D [cu 2 + (fc 2 D) 2 ] -1 \, and the current-current correlator, 
Ajiji(u>,k) — 2uj 2 S tl D[uj 2 + (fc 2 £>) 2 ] _1 x. As required, this current-current correlator is consistent with 
the Kubo formula of footnote 11. 
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Analogous results for the real-time correlations of energy and momentum densities are 

A ££ {t,k) = e -5 fc2 7*l*l cos(\k\v s t) C , (33a) 

iA ni£ (t,k) = F^e"5 fc Vl*l S m{\k\v a t) C, (33b) 

iA mS {t,k) = (W/v^e-^W sm(\k\v s t) <|P 2 )/V, (33c) 

A^{t,k) = UP' - Pfe) e-*^W + k'P e-v*' 1 ''® cos(|fc|u s t)l {\P 2 ) /V . (33d) 

where (£, fc) = f d 3 x e~ tk x {7r*(t, a?), 7r J '(0, 0)}), etc, and C is the mean square 
fluctuation in energy per unit volume defined in Eq. (17). The mean-square fluctuation 
in total spatial momentum may be deduced from the equilibrium relation (11) between 
momentum density and flow velocity, which implies 

(pi pi) i Qt? 

VSft7 = * r ' u ' (34) 

or (|P 2 )/V = wT. Hence the mean square fluctuations in momentum and energy are 
related by the velocity of sound, 

hfLl = Tw = T 2 ^- = v 2 s T 2 ^-=v 2 C. (35) 
V dT s dT y J 

E. Long-time tails 

Consider a spatial charge flux integrated over all space, 

J a {t) = J d 3 xj a (t,x). (36) 

If non-linear terms in the constitutive relation for j a are ignored, this zero wavenumber 
charge flux has no coupling to hydrodynamic variables due to the gradient in the linear 
—DVii part of its constitutive relation. Hence, a linear response analysis predicts that the 
zero wavenumber flux should relax on a short, microscopic time scale. 20 This conclusion is 
wrong, because the charge flux does couple to hydrodynamic degrees of freedom, even at zero 
wavenumber, through the non-linear niz/w term in its constitutive relation (8). Although 
this term is quadratic in deviations away from the chosen equilibrium state, because it 
contains no spatial gradients it will always dominate over the single gradient —DVn term 
of linear response for sufficiently long wavelengths. 

To evaluate the zero-wavenumber (J*(t) J^(0)) correlator correctly for large times, one 
needs only to express these currents in terms of the product of charge and momentum 
density fluctuations, and then insert the results (30) and (33d) for the spectrum of these 



20 This may be seen explicitly from the linear response form of the current-current correlator A tl (uj, k) shown 
in footnote 19, whose k — > limit is frequency independent, lim^o A ll (u>, k) = 25 ll D\. 
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fluctuations: 



= Hpi [ AJL [| e - fc2 (^)l*l x + I e - fc2 ( D +^)l*lcos(|fc|^) X " 
w J {2ix) 6 L d 6 

; a + (exponent., decay) . 




rf 3 x (|{n a (i,xy(*,:r), n 6 (0, 0)^(0,0)}) 
rf 3 x <|K(i, x), n 6 (0, 0)}) x), 7T J '(0, 0)}) 



- aft 



(37) 



Factorizing the first correlation function into a product of two-point correlators is justified 
because, in the small-fe hydrodynamic regime, the distribution of fluctuations is arbitrarily 
close to Gaussian. 21 After inserting the explicit forms of C narib (t,k) and C w i n j(t,k) in the 
third step, an angular average over the direction of k was performed. In the final step, the 
term involving longitudinal momentum fluctuations does not contribute to the t -3 / 2 long- 
time tail because of the oscillating cosine; this term only gives rise to an exponentially falling 
0[exp{ — ||t|^s/(-D+|7s)}] contribution. The above treatment is essentially perturbative 
theory for small fluctuations in the long-lived degrees of freedom; upon factorization, the 
two-point correlators of conserved densities are taken to evolve according to linear response. 

A completely analogous long-time tail must also be present in the zero wavenumber stress- 
stress correlator due to the coupling of the stress tensor to products of energy, momentum, 
and charge density fluctuations. Inserting the constitutive relation (21), and evaluating all 
the terms which can contribute to a long-time tail is straightforward, but slightly tedious. 
This calculation is summarized in Appendix A. Here we will specialize to the simpler case 
of a scale invariant theory, for which 



The essential point is that hydrodynamic variables represent microscopic fields averaged over spatial 
volumes which are large compared to the "correlation volume" £ 3 (where £ is the static correlation length) . 
Hence, a hydrodynamic fluctuation with wavelength A 3> £ may be regarded as averaging over roughly 
(A/£) 3 essentially independent correlation volumes. By the central limit theorem, the distribution of such 
an average will approach a Gaussian distribution as the number of independent samples becomes large. 
Residual non-Gaussian correlations in fluctuations with wavenumbers of order k will generate corrections 
to the result (37) which arc suppressed by an additional factor of |fc| 3 £ 3 , leading to sub-dominant 0(t~ 3 ) 
corrections to the final result. 



where E^ n 




\ 5 l i 5 mn is a projector onto traceless symmetric tensors. 



V + \ 5e(t, x)} + - W± vr m (t, x) n n (t, x) ) , (38) 
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Therefore the zero-wavenumber stress-stress correlator in the hydrodynamic regime is 
V-^HT^T^O)}) = I d 3 x j Hf q {\ {vr m (t, x) n n (t, x), tt^O, 0) tt«(0, 0)}> 

+ I^ < 5«(I{ e (*, a .), e (0,0)})J. (39) 

The first term may be evaluated using the same logic as above. Retaining only pieces which 
contribute to the connected correlator [i.e. (|{T lJ (t), T kl (0)}) — (T y '(i))(T fcZ (0))], one finds 

^ / d 3 xHH n H k p l q (l{n m (t,x)n n (t,x), ^(0,0)^(0,0)})^ 

= J^flJJ j d 3 x (l{7r m (t,x),^(0,0)})(l{7r n (t,x),^(0,0)}) 
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(i)3/2+ 1ra^' (40) 

Once again, contributions which fall exponentially with time have been dropped in the last 
step. The final term in Eq. (39) gives a time-independent contribution, 

\5 lj 5 kl J d 3 x (±{e(t,a:), £(0,0)}) = \ 5 lj 5 kl A ££ {t,k=0) = ±eT5 ij 5 kl . (41) 

This reflects the direct coupling between stress and fluctuations in energy density. At zero 
wavenumber, these become fluctuations in the total energy (in the equilibrium grand canon- 
ical ensemble), which are strictly conserved in time. Hence, one finds that the leading long- 
time behavior of the connected zero-wavenumber stress-stress correlator in a scale invariant 
theory is 22 



+ 15 



U + 1 \ b° + 2° 3° ) ( 8 7r 7?? |t|) 3 / 2 • 

(42) 



Adding any higher-order terms to the stress tensor constitutive relation (23) will generate 
sub-dominant power- law tails in (42) proportional to |t|~ 5 / 2 , |t|~ 7 / 2 , etc. This is because such 
higher-order terms will contain either additional spatial gradients or additional powers of 
fluctuations away from equilibrium. In either case, one finds contributions which scale with 
higher powers of A; as k — > 0, implying additional powers of l/\t\. 



Once again, this result is for three spatial dimensions. In d > 3 spatial dimensions, the same analysis 
leads to a |i|~ d / 2 long-time tail. In two or fewer spatial dimensions, hydrodynamic fluctuations are 
sufficiently infrared-singular that the oj — > limits in the Kubo formulas in footnotes 11 and 14 fail to 
exist. Instead, one finds that transport coefficients exhibit non-trivial scale dependence on arbitrarily long 
length scales [18]. 
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F. Low frequency behavior 



The long-time tails present in the current-current correlator (37) or the stress-stress 
correlator (42) automatically imply that the Fourier transforms of these correlation functions 
cannot be analytic at zero frequency. 23 A little analysis shows that if G(t) ~ a /\t\ 3/2 at large 
times then 24 

G(u) = J dt e iuJt G(t) V^a\uo\ 1/2 + (analytic terms) , (43) 

as uj — > 0. Consequently, in three spatial dimensions, the Fourier transforms of the real-time 
zero-wavenumber correlators (§{•/„(£), J 3 b (0)}) and (|{T* J (t), T M (0)}) must have square 
root branch points at zero frequency whose coefficients are related via Eq. (43) to the 
amplitudes of their long-time tails. 

In thermal equilibrium, the spectral density of any two operators A and B, 

p AB (u) = J dte™ ([A(t),B(0)\) , (44) 

is directly related to the Fourier transform of the corresponding symmetrized correlator 

Pab (uj) = 2tanh(/^/2) J dt e iu)t (\{A(t),B(0)}) . (45) 

Consequently, the current-current or stress-stress spectral densities, divided by f3u, must 
have the same non-analyticity at uj = as do the Fourier transformed real-time correlators. 
For example, in a scale-invariant theory one has 

- / dx°d 3 xe^° ([T y (x),T H (0)]> = (I + 

x [sW + s^M^ + OtM 3 / 2 )], (46) 

with 



s (o) = 4?7> s (i) 



60tt 



G. Large iV c scaling 

The stress-stress correlator will, of course, have short time transient contributions in 
addition to its long time power-law tail. In large N c gauge theories (or large iV matrix 



This follows from the standard argument showing that if an integrable function G(w) is analytic in a strip 
of width 25 surrounding the real axis, then its Fourier transform G(t) falls off exponentially as 0(e^ 5 '*'), 
as may be seen by shifting the contour of integration. 

More generally, if G(t) ~ a|t|~ d//2 , then the leading non-analyticity in G(u>) for small frequency is G{w) ~ 
a [i(d-i)(d+i)/4 V2 7r/r(d/2)] \u\( d - 2 ^ 2 for d odd, G(u) - a [7r(-l) d / 4 /r(d/2)] |w|< d - 2 )/ 2 for d = An, and 
G(lu) - a [(-l)(<i+2)/4 2/r(d/2)] l^l^- 2 )/ 2 In |w| for d = 4n + 2. 
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models) if the 't Hooft coupling g 2 N c is held fixed as N c becomes large, then the short time 
contribution to the (connected) stress-stress correlator will scale as iV 2 , just because there 
are 0(iV 2 ) degrees of freedom. In contrast, from Eq. (42) one sees that the coefficient of the 
£-3/2 j on g tj me t a jj w ill be 0(1) as N c — > oo if the shear viscosity divided by the enthalpy, 
•y^ = rj/w, has a finite, non-zero large N c limit. This ratio is effectively a microscopic 
relaxation time, the mean free time for large angle scattering of elementary excitations 
in the system. In weakly coupled high temperature gauge theories, this ratio does have 
a finite large N c limit [19]. Presumably this is generally true, at least in theories like 
Af—A supersymmetric Yang-Mills theory which are strongly believed not to have any phase 
transition separating their weak and strong coupling regimes. Assuming this is the case, the 
cross-over time (beyond which the long-time tail dominates over the short time transients) 
will become arbitrarily large as N c — > oo. 25 For any fixed time t, the t~ 3//2 tail will not be 
visible in the leading large N c behavior of the stress-stress correlator, but will instead appear 
as a sub-leading 1/iV 2 relative correction: 

(T y (t)T„(0)> ~ 0(N 2 C ) e-*' T + 0(1) r 8 ' 2 . 

Exactly the same result is true of the long-time tail in current-current correlators. Diffusion 
constants are 0(1) as N c — > oo (at least in weakly coupled hot gauge theories [19]) and 
therefore for fixed time t, the long time tail in the current correlator (37) also scales as 
0(l/iV 2 ) relative to the short time transients (which scale the same as the susceptibility x)- 
In both cases, the result that long-time tails become sub-leading effects at large N c 
may be regarded as a consequence of the fact that equilibrium velocity fluctuations are 
"anomalously" small when N c — > oo. As noted earlier in section II B, the flow velocity equals 
the momentum density divided by the enthalpy. Therefore, the mean square fluctuations in 
average flow velocity are directly related to the mean square fluctuations in total momentum. 
Using the previous result (34), one has 

/ pi pj\ Xij t 

W = W=«v = 0(J ^ (48) 

as N c — > oo for fixed volume. More generally, this means that the fluctuations in flow 
velocity averaged over a patch of fluid of linear size L have a 1/N C suppression (on top of the 
expected L~ 3 / 2 behavior). Consequently, in the spatial current density j = —D Vn+n v, the 
non-linear nv contribution is 1/N C suppressed relative to the linear —D'Vn term. Exactly 
the same conclusion holds for the terms quadratic in flow velocity in the stress tensor (23) 
as compared to the linear terms involving the gradient of the velocity. 



Completely analogous nonuniformities between the large distance and large N limits are well known in 
other contexts. For example, adjoint representation Wilson loops in SU(A r c ) Yang-Mills theory have 
a cross-over from area- law to perimeter- law behavior on a distance scale which diverges as iV c — ► oo. 
And various SU(N) symmetric two-dimensional models have correlators which fall with distance like 
\x\-V N [20]. 
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III. SUPERSYMMETRIC HYDRODYNAMICS 



We now wish to generalize the preceding discussion of hydrodynamic fluctuations to the 
case of supersymmetric theories. For ease of presentation, in this section we consider theories 
with J\f—1 supersymmetry in four dimensions. Such theories possess a set of conserved 
super currents, which we will denote by S£(x) and S£(x) , which are Lorentz vector-spinors. 
The corresponding conserved supercharges satisfy the supersymmetry algebra, 26 

{Q f ,,Qi} = 2afaP li . (49) 

The theory may also have a chiral U(l) symmetry (i?-symmetry) , whose charge R does not 
commute with the supercharges, 

[Qa,R}=Qa, [Qa,R] = -Qa. (50) 

as well as various additional internal symmetries whose charges N a do commute with the 
supercharges. For definiteness, we will assume that the theory does have i?-symmetry, and 
that this symmetry is not spontaneously broken. 27 

Since the supercurrents are conserved, one might expect that fluctuations in supercharge 
densities will be unable to relax locally, just as ordinary current conservation directly leads 
to diffusive relaxation of bosonic conserved charge densities. And slow relaxation of super- 
charge densities may in turn generate long-time tails in bosonic observables which can couple 
quadratically to supercharge densities. In other words, even though supercharge densities 
are fermionic operators with vanishing expectation value in any conventional statistical en- 
semble, the dynamics of supercharge fluctuations may produce distinctive effects in ordinary 
bosonic observables. We will see that these expectations are largely correct. 



A. Supercharge fluctuations 

The hydrodynamic nature of supercharge densities is reflected in the long-time behavior 
of the corresponding correlation function, 

CZ{x) = (S&x) S£(0)) . (51) 



We follow the notations of Wess and Bagger [21] for spinors. In particular, undotted early Greek indices 
a, /?,... label (§, 0) two-component Weyl spinors, while corresponding dotted indices label (0, 5) conjugate 
Weyl spinors. These indices are raised or lowered using the antisymmetric tensor e 2 i = e 12 = e 2 i = 
e i2 = 1, so that £, a Q a = -£, a Q a and £ & Q & = -£aQ & - Extended Pauli matrices are 0* ■ = (-l,cr) 
and <j^ a P = (—l,—o"), and the spinor representation generators of the Lorentz group are (<J^ V )£ = 
\{<a^ <^ & ) and = ~ °™ a ^J- 

In 7V=1 supersymmetry, conservation of the U(l) i?-symmetry current requires scale invariancc. Non- 
scale invariant theories will have anomalies in the i?-symmetry current, although it is often then possible 
to combine the anomalous i?-symmctry with an anomalous chiral symmetry to yield a non-anomalous 
i?-symmetry. If there is no non-anomalous i?-symmetry, then the i?-charge density will not be a hydro- 
dynamic degree of freedom For such theories, just ignore all mention of the i?-symmetry current in the 
following discussion. 
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Finite-temperature thermal ensembles are not invariant under supersymmetry transforma- 
tions. 28 This is reflected in the fact that at finite temperature bosons and fermions obey 
different statistics (Bose-Einstein or Fermi-Dirac), as well as in the fact that bosonic and 
fermionic fields obey different boundary conditions (periodic versus anti-periodic) in Eu- 
clidean functional integrals representing equilibrium thermal systems. But it is important 
to understand that the non-invariance of thermal ensembles under supersymmetry does not 
imply that the supercurrent fails to be conserved. The conservation of the supercurrent, 

dMx) = (52) 

is an operator identity, valid in any physical state. Since thermal averages are just linear 
combinations of expectation values in physical states, the conservation of the supercurrent 
implies that the thermal supercurrent correlation function satisfies the Ward identity 29 

d,C^(x) = 0. (53) 

Similarly, taking thermal expectations of both sides of the supersymmetry algebra (49) shows 
that mean square fluctuations in supercharge are directly related to the total energy of the 
thermal ensemble. A spatial Fourier transform of the Ward identity (53) gives 

d Cg(fc, t) = -ik 3 Cg(fc, t) = 0{k) , (54) 

which directly shows that long wavelength fluctuations in supercharge density must relax 
arbitrarily slowly. 



B. Super-thermodynamics 

Having realized that supercharge densities are hydrodynamic degrees of freedom in su- 
persymmetric theories, we would like to construct the appropriate hydrodynamic equations 



Formally, the transformation of a statistical density matrix p under any infinitesimal symmetry trans- 
formation is defined by the condition that Tv(Sp O) = Tr(p SO), where O is any observable. This is a 
generalization of the relation between Schrodinger and Heisenberg pictures, and expresses the equivalence 
between active and passive views of symmetry transformations. If Q is the generator of the transformation, 
then the infinitesimal transformation of the observable is SO = [Q,0]. For normal (bosonic) symmetries, 
this of course implies that 5p = —[G,p], so invariance under the symmetry means the density matrix 
p commutes with the symmetry generator Q. But for an infinitesimal supersymmetry transformation, 
G = (, a Qa w hh £ an external Grassmann spinor which anticommutes with all fermionic operators. In 
this case, one has Tr(p [£Q, O]) = Tr(£/> [Q, O]) = Tr([£p, Q] O) = Tr({p, £Q} O), so the variation of the 
density matrix is the anticommutator with the supersymmetry generator, 5p = {p,^Q} = £,{p,Q}, not 
the commutator. [We assume throughout that p and O are bosonic operators — possibly constructed by 
multiplying odd numbers of fermionic operators by independent external Grassmann spinors.] Therefore, 
even though the equilibrium canonical ensemble p = Z~ x e~$ H commutes with the supercharge Q, this 
does not mean that the canonical ensemble is invariant under supersymmetry. 

This assumes one suitably defines the short-distance subtraction in the renormalization of the product of 
supercurrents to remove what would otherwise be contact terms on the right hand side. One can always 
do so, provided the supersymmetric theory under discussion exists. 
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in complete analogy with the earlier non-supersymmetric treatment. To do so, we must 
construct the appropriate constitutive relations for the spatial part of the supercurrent (and 
the /^-symmetry current), and also understand possible new terms involving products of 
superchange densities which might appear in the constitutive relations for the stress tensor 
and other ordinary symmetry currents. 

As discussed previously, non-derivative terms in constitutive relations are completely de- 
termined by equilibrium thermodynamic derivatives of the free energy. This is true provided 
one has introduced chemical potentials conjugate to all conserved densities of interest. To 
do this for a supersymmetric theory, we must first understand the generalization of ordinary 
thermodynamics to 'equilibrium' ensembles in which there is a non-zero expectation value 
of the supercharge, produced by turning on a super-chemical potential. 

To construct an equilibrium state with non-zero supercharge density, consider generalizing 
the usual statistical operator (3) by adding supercharges multiplied by Grassmann-valued 
super-chemical potentials, 30 

p s = Z- 1 exp [0 (u v P v + p a Q a + p a Q«)} . (55) 

To simplify expressions, we are temporarily suppressing any chemical potentials for bosonic 
conserved charges. Note that the addition of spinorial chemical potentials coupled to the 
supercharges Q a and Q a maintains translation invariance of the equilibrium state but ex- 
plicitly breaks rotation invariance. 

In the ensemble described by p s , the supercharges Q a and Q a have non-zero expectation 
values proportional to p and p, respectively. More precisely, to first order in p and p, 

= -P(p & QaCQa)/V 

= (5C{ea°-wa-u)^p^ (56) 

where 0{u 2 ) corrections have been dropped in the last equality. Here £ is an independent 
Grassmann spinor. The intermediate expectations in (56) denote expectation values in the 



Though fcrmionic chemical potentials may seem peculiar, one can regard them as a non-dynamical grav- 
itino field, to which the theory is coupled. It should be stressed that these generalized density operators 
do not have the usual statistical interpretation. Their diagonal matrix elements, although commuting, 
are not real numbers, and therefore cannot be viewed as statistical probabilities. Nevertheless, these 
generalized density operators may formally be regarded as equilibrium ensembles, in which generalized 
thermodynamic relations hold. The logical alternative of introducing real-valued chemical potentials for 
supercharges is formally inconvenient and also fails to yield any satisfactory physical interpretation. In 
particular, a density matrix containing a real-valued chemical potential for supercharges, when applied 
to the union of two disjoint physical systems, does not factorize, even if the two systems are causally 
disconnected [22]. 
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usual canonical ensemble (with 4- velocity u), and trace cyclicity and the supersymmentry 
algebra (49) have been used. Expanding to higher order in p and p, would lead to corrections 
of order pp 2 which will not be relevant for our purposes. Similarly, 



up to corrections of order p 2 p [and 0{u 2 )}. 

Naively, one would expect the generalized partition function 



(57) 



Z(0, u, p) = Tr exp [{3 (u v P v + p a Q a + /2 d Q*)] (58) 

to be a generating function for expectation values of supercharges and products of super- 
charges. However, this is not true. Although the two terms we have added to the exponential 
defining p s both commute with the total momentum P u , they do not commute with each 
other. As a result, the derivative (d/dp) In Z{f3, u, At, At) does not equal P{Q). In fact, the 
generalized partition function is completely independent of the super-chemical potentials, 
and simply equals the usual (physical) partition function, 

Z(P, u, At, At) = Z((3, u, 0, 0) = Tr e pUvPV . (59) 

To see this, first note that 

Tr [p a e^ pv ■ ■ ■ ] = p a Tr [(-1) F e^ p " • • • ] , (60) 

and therefore 



9 -Tr[e^...]=Tr[(-l)^ f 



■]■ 



(61) 



dp a 1 dp a 

Here (— 1) F = e 2ntJz is the operator which multiplies all bosonic states by +1 and all 
fermionic states by —1. This operator commutes with P u and anticommutes with the 
supercharges. Next, factorize the exponentials in p s using the Baker-Campbell-Hausdorff 
formula, 



(62) 



Consequently, 



1 d 



(3 dp~ 



Z(P,u,p,p) = Ilr [(-1) f ^ ^°«« e^^-^"^ 4 ) 



= Tr 
= Tr 
= Tr 



e to°Q a e ^Q* JPW-W <* d P d ) (_g 7 _ p<%. jf Pv ) 
(-If (-Q 7 + W & Ji & P v ) e^ aQ « e PMu»-^K^) 1 . ( 63 ) 



The penultimate step uses the anticommutation of Q with (— 1) F and trace cyclicity to move 
Q from the left to the right end of the trace. The final step uses the supersymmetry algebra 
to commute Q through the exponential of Q: 



^7 j c- 



= PjP{Qv Q-y} e M& = 2(3 P v . 



(64) 
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Since the last line in (63) is minus the second line, this shows that ( y d/dfi' y )Z(f3, u, //, p) = 0. 
One may show that (d/dp,^)Z(/3, u, //, p) = in a completely similar fashion. 

In other words, although adding fermionic chemical potentials to the statistical opera- 
tor induces non-zero expectations values for the supercharge densities, it does not change 
the partition function at all. This means that the "thermodynamic pressure", defined as 
(/9V) _1 In Z, is not affected by the presence of non-zero equilibrium supercharge densities. In 
fact, this is required if the thermodynamic pressure is to agree with the microscopic pressure, 
defined in terms of the expectation value of the stress tensor, V = |Tr (p s Tn). To see that 
the microscopic pressure is also /z-independent, one may repeat the steps shown in Eq. (63) 
when there is an additional insertion of the stress-energy tensor. This yields 

^^W = H(-V F [Q°>T»])^. (65) 

But the supersymmetry transformation of the stress-energy tensor involves spacetime deriva- 
tives of the supercurrent (see, e.g., [23]), 

lQ a , th = -| { ko/ d p s» + (o / dpSf] > (66) 

whose expectation (65) vanishes in the translationally invariant equilibrium ensemble. 

None of the above results are affected by the presence of non-zero chemical potentials for 
ordinary conserved charges which commute the the supercharges. But if a non-zero chemical 
potential /j, R for the .R-charge is present, then the generalized partition function does acquire 
dependence on the fermionic chemical potentials. The lowest such term is proportional to 
[ififi R ; this term is directly related to the equilibrium form of the -R-charge current shown 
below. 

In a similar fashion, one may evaluate the generalized equilibrium expectation values of 
the bosonic symmetry currents and the supercurrents. Ordinary bosonic conserved currents 
commute with the supercharges, and hence are unaffected by super-chemical potentials. In 
the presence of non-zero ordinary chemical potentials, 

<^U = ^ (67) 

as usual, where n a is the equilibrium density of the conserved charge N a . For the supercur- 
rents one finds, 

(CS x a )^ = -^nr^^+OW, (68a) 

<^«a>^ = ^n^^a+O^), (68b) 

where T^ v = w u^u v + V g^ v is the usual equilibrium form of the stress-energy tensor. These 
results reflect the supersymmetry transformations 

{Qp, S*} = 2T\ + (space-time gradients) , (69a) 
{Q^, S^} = 2T\ + (space-time gradients) , (69b) 

and {Qp, S*}, {Qp,S£} being pure space-time gradients. (The expectation value of the 
space-time gradient terms vanish in all translationally invariant equilibrium states.) For the 
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A = components of these commutators, the unspecified space-time gradients are strictly 
spatial gradients which vanish on spatial integration, so that the anti-commutation relations 
(69) are consistent with the supersymmetry algebra (49). For the i?-symmetry current, one 
finds 

< J R \_ = n R u x -\p f\ (^ f) + 0(/iV) , (70) 

where n R is the equilibrium i?-charge density in the absence of super-chemical potentials. 
The term quadratic in super-chemical potentials arises from the non-trivial commutators of 
the supercharges with the .R-current, 

[Qa, Jr\ — S£ + (space-time gradients) , (71a) 
[Qa, Jr\ = —S^ + (space-time gradients) . (71b) 

Note that the relation (70) implies that non-zero super-chemical potentials induce a non-zero 
i?-charge density even when the -R-charge chemical potential vanishes. 



C. Super symmetric constitutive relations 



We can now construct the required constitutive relations for the stress tensor and super- 
current. As before, we will need to include all possible terms which are linear in fluctuations 
away from a reference equilibrium state (chosen to have vanishing charge, supercharge, and 
spatial momentum densities) involving one spatial derivative, as well as terms without deriva- 
tives which are linear or quadratic in fluctuations. Neglected terms involving either more 
derivatives or more fluctuations will not affect the leading long-time tails in stress-stress or 
current-current correlators. 

We will use p a and p a to denote the (Grassmann valued) supercharge densities in the 
non-equilibrium ensemble of interest. In equilibrium, the relations (68) imply that the 
supercharge densities are related to the super-chemical potentials via 



Pc 



w a 



Pa 



-/?//( 



ea 



w a 



u ) a- 

I pa 



(72) 



neglecting 0(u 2 ) corrections. For the reasons discussed in the previous subsection, to 
quadratic order in deviations away from the reference equilibrium state, the pressure con- 
tains no coupling to supercharge densities. Consequently, the constitutive relation for the 
stress is identical to the previous result (21), except for the inclusion of i?-charge density 
fluctuations along with other bosonic charge fluctuations, 31 



V + v 2 5e + if (5e) 2 + \ E ab n a n b + \ S R n^ 2 - v 2 — 



w 



+ 



71 71 
W 



- 7 C 5 ij V • 7T - 7t; (V V + VV - p ij V • tt) . 



(73) 



31 In writing the constitutive relation (73), we have assumed, for simplicity, that the global bosonic symmetry 
currents j a are all vectors, rather than pseudo- vectors, and that parity is a symmetry of the theory. This 
prevents a S aR n a n R cross term from appearing in the constitutive relation (73) for the stress, and likewise 
ensures that the i?-charge density does not mix with other global charge densities in the constitutive 
relations (74) and (75) for the currents. 
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The constitutive relations for bosonic internal symmetry currents are unchanged from 
the non-supersymmetric case, 

11 7T 

j a = ~D ab Vn b + ^- , (74) 

w 

while the new constitutive relation for the i?-symmetry current is 

71 IT 

j R = - A Vn R + ^- - — p„ Pp . (75) 
w 2e z 

The quadratic coupling to supercharge densities follows directly from the (generalized) equi- 
librium form (70) for the expectation of the i?-current, combined with the relation (72). 

The remaining constitutive relations are those of the supercurrents. The appropriate form 
is more complicated than for bosonic currents because the supercharge density is a spinor. 
There are two independent rotationally covariant forms in which a spatial gradient can be 
applied to a spinor p a to form a vector-spinor, namely V l p a and (cr l *)(f V,- pp. Consequently, 
the supercurrent constitutive relations have the form 

Si, = -D&Pa-D^fVjPfi 

Si = -D s V% - D.Vj p s (a j % 

Here D s , D a are two new diffusion constants, which multiply the two structures allowed by 
rotation invariance. 32 The remaining non-derivative terms are dictated by the equilibrium 
results (68), when expanded to first order in flow velocity. The third term, which is linear in 
supercharge density, reflects the fact that a constant supercharge density is not rotationally 
invariant, and hence induces a non-zero constant spatial supercurrent. The convective terms 
are more complicated than just pa^/w, because of the presence of this linear non-derivative 
term. If the supersymmetric theory in question is also scale invariant, then a^S^ must 
vanish, as it is in the same supersymmetry multiplet as the trace of the energy-momentum 
tensor. Applying this constraint to the constitutive relation (76), one finds that the two 
diffusion constants must coincide, D a — D s , in scale-invariant supersymmetric theories. 



- | V P, + \ Pa ^ + j 2 ttV* V pp , (76a) 

- | Pp (SV)^ + \p^ 1 + | (v k * l A ■ (76b) 



D. Relaxation of supercharge fluctuations 

The constitutive relation (76), combined with conservation of the supercurrent, allows one 
to determine how long-wavelength fluctuations in supercharge density relax. Inserting the 
constitutive relation into the continuity equation d t p a = — V-<S a , neglecting the non-linear 
convective terms, and Fourier transforming in space yields 

[(d t + D 8 k 2 ) Si - ic s (k ■ a a°)f] p p (t, k) = 0, (77) 



Assuming parity invariance of the theory ensures that the diffusion constants for the left and right handed 
supercurrents are the same. The Kubo formula for the supercharge diffusion constant D s is e D s = 
- lim^o j2 J d 4 xe iux °a ^ Cf fj {x). 
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where c s = V/e. The solution in the small-fc limit shows weakly damped, propagating 
behavior, 



Pa(t, k) = e 



-D s k 2 t 



5 P COS (\k\c s t) +i(k- era )/ sin(|fc|c s t) P/3 (0, k) 



(78) 



This is a propagating collective excitation which is distinct from normal sound. In some 
respects, it is analogous to second sound in a superfluid. Note that the velocity c s = V/e 
of these "supersound waves" is always less than the ordinary speed of sound t> s = a/ dV / de. 
In scale invariant theories, c s = 1/3 while v s = 1/V3. Collective excitations in supercharge 
density were discussed earlier in Ref. [24], from a somewhat different perspective. 

Just as in the previous non-supersymmetric analysis, the linearized non-equilibrium relax- 
ation (78) may be converted into a statement about the time-dependence of the equilibrium 
supercharge density correlator. Let 



C(f,fc)^ = J d 3 xe- ik - x (S a (t,x)S%0,0)), 



(79) 



where the expectation value is in the reference equilibrium state with vanishing chemical 
potentials. The small k limit of the equal time correlator is fixed by the supersymmetry 
algebra, 33 

C7(0,fc=0)^ = -ea^. (80) 

Consequently, the linear response result (78) is equivalent to the statement that, in the small 
k limit, the equilibrium supercharge correlator for large times has the form 



C (^ k ) a p = -ee 



-D s k 2 \t\ 



a° a cos(\k\c s t) + i(k ■ cr) a ^ sin(|fc|c s t) 



(81) 



One finds the same result for the small k behavior of the supercharge correlator with the 
opposite ordering, 



C(t,k) ap ee J d 3 xe-* kx (S°(t,x)S a (0,0)) 



a° a/ j cos(|fc|c s t) + i(k ■ <r) a p sin(|fc|c s t) 



(82) 



E. Supersymmetric long-time tails 

As in the non-supersymmetric case, second order terms in the constitutive relations are 
sufficient to understand the leading long-time tails in k=0 correlators of all conserved cur- 
rents. Proceeding in exactly the same way as in Section II E, one sees from the absence of any 
pp term in the stress tensor constitutive relation (73) that supercharge density fluctuations 
have no effect on the leading long-time tail in the stress-stress correlator. Consequently, the 
stress-stress correlator has the same form as in the non-supersymmetric case [i.e., Eq. (42) for 



33 In our conventions, cr° = — 1, so C(t=0, k=0) is strictly positive, as the equal-time correlation function of 
any operator with its Hermitian conjugate must be. 
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scale-invariant theories, or Eq. (All) for the general case] provided the i?-charge is included 
as one of the global bosonic charges. Similarly, the leading long-time tails (37) in correlators 
of ordinary bosonic symmetry currents are unaffected by the presence of supercharge density 
fluctuations. 

In contrast, the presence of a quadratic supercharge density term in the constitutive 
relation (75) for the i?-symmetry current means that the long time tail in the i?-symmetry 
current correlator is affected by supercharge fluctuations. One finds that the leading large 
time behavior is 

/*» m^UM}) ~ £ { pgg^p + • <*> 

If the same analysis is applied to the spatial supercurrent correlation function, one finds 
that the second order irp terms in the constitutive relation (76) do not generate long-time 
tails. Because the speed of ordinary sound (v s ) and supersound (c s ) differ, fluctuations 
(with the same wavevector) in momentum density and supercharge density cannot remain in 
phase. This is analogous to the fact that only transverse momentum density fluctuations, not 
longitudinal sound waves, contribute to the long-time tail (37) of ordinary current-current 
correlators. So the supercurrent correlator {S^it) S^(0)) , at k=0, should decay exponentially 
in time. 34 

The differing long-time behaviors of the stress tensor and supercurrent correlators does 
not conflict with the fact that the stress-energy tensor and supercurrents belong to the 
same supersymmetry multiplet. The supersymmetry algebra implies simple relations be- 
tween (— 1) F inserted thermal correlation functions of conserved currents, but not between 
conventional thermal correlation functions, for the reasons noted in footnote 28. 



IV. LONG-TIME TAILS IN M = 4 SUPERSYMMETRIC YANG-MILLS THEORY 

As a specific application, the above discussion (trivially generalized) may be applied to 
jV=4 supersymmetric SU(N C ) Yang-Mills theory, at the origin of the moduli space and 
at non-zero temperature. 35 Once again, hydrodynamic degrees of freedom are densities of 
conserved charges. In Af—4 supersymmetric Yang-Mills theory, the supersymmetry algebra 36 
is generated by the charges of the dilation current jp = x v T^ u , the special conformal currents 
K^^ = 2x\x u T^ v — x 2 T' I x , and the special supersymmetry currents S^ a = x u^aa^T anc ^ 
^Aa = x v ®6ia ^a* > m addition to the stress-energy tensor T^ u and supercurrents S^ a , S^. 
The index A = 1,...,4 is the extended supersymmetry label. There are also 15 bosonic 
i?-symmetry currents, whose charges generate a non-anomalous global SU(4:) R symmetry. 



Because c s < v s , one finds that higher-order pir n terms in the supercurrent constitutive relation also do 
not generate (higher-order) power-law tails. 

The vacuum degeneracy associated with a non-trivial moduli space is lifted at non-zero temperature. For 
sufficiently high temperatures, the only stable equilibrium states lie at the origin of moduli space. 
See, for example, Ref. [23], for an introduction to the 7V=4 algebra. 
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The currents j^, i^ M A , S^ a , and S^a are all constructed from the energy- momentum tensor 
and ordinary supercurrents. Therefore, the time evolution of the charge densities of these 
currents is completely determined by the time evolution of the energy-momentum tensor 
and the supercurrent. In other words, j° D) K° x , S Aa , and <S A6l are not new hydrodynamic 
degrees of freedom, independent from T 0l/ , S Aa , and S Aa . No Abelian magnetic fields appear 
as hydrodynamic degrees of freedom, since the gauge group of the theory is SU(N C ), not 
U (N c ). Moreover, holographic arguments imply that the theory at non-zero temperature has 
a mass gap [5], in accord with the assumption of finite correlation length made in Section II. 

The discussion of supercharge fluctuations in the previous section generalizes easily to the 
case of Af—4: supersymmetry. In particular, the constitutive relation for the stress tensor 
has the same form (73), except for the addition of a fundamental SU(4:)r index on the 
supercharges p Aa and an adjoint index on the i?-charge densities n Ra (i.e., a = 1, 15). 
The theory is scale- invariant, and therefore e = 3V, and the bulk viscosity vanishes, as 
does the coefficient E R . Hence, fluctuations in the i?-charge densities do not contribute to 
constitutive relation for stress to second order, which therefore must have the same form 
(23) as in a non-supersymmetric scale-invariant theory. Consequently, the long-time tail for 
the stress tensor correlator (at fc=0) is given by the previous result (42), and the equivalent 
small frequency form of the spectral density is given by Eq. (46). 

The constitutive relation for supercurrent (76) stays unchanged to second order (except 
for the addition of the index A, and the specialization to D s = D a and P/e = 1/3). No long- 
time tail is present in the supercurrent correlator, for the reasons discussed in the previous 
section. Hence, the corresponding spectral density is analytic at zero frequency. 

The constitutive relations for the i?-symmetry currents become 

j Ra = -£> R Vrfe a + - PAa {ta) AB d** P B P , (84) 

w 6e 

where t a are the generators of SU(A), normalized so that tr(t a t b ) = \8 ah . The resulting 
long-time tails in the i?-symmetry current correlators (at k=0) are 

j** m.(M(o,o) } ) ~ [W ,g: wl 3 /2 + (85) 

The equivalent small-frequency expansion of the spectral density is given by 

i J dtd 3 x e™(\£ a (t,x),£ h (0M) = W (c (0) + c (1) VF\) +0(M 3 / 2 ) , (86) 

where 

c (o) 

In writing expressions (84) and (85) we have used the fact that in a thermal state without 
chemical potentials for the SU (A) R charges, the matrix of diffusion constants, as well as the 
R charge susceptibility, must be invariants of the group: (£k) a & = D R 5 ab , (xnjab = Xk<U- 



2MiXr ion \ 

-^r- , (87a) 
-g-M^ )Y *» D?». (87b) 
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The amplitudes of the long-time tails for the stress tensor correlator (42), and the R- 
current correlator (85), depend on the values of the equilibrium energy density e, the R- 
charge susceptibility Xr? the shear viscosity r], and i?-charge diffusion constant D R . These 
parameters must be treated as input from short-distance physics. In a weakly coupled theory, 
these quantities can be computed in perturbation theory, but no field-theoretic calculation 
is available when the coupling constant is large. However, the AdS/CFT correspondence 
predicts specific values for these parameters in J\f—A supersymmetric Yang-Mills theory in 
the limit of large iV c and large 't Hooft coupling. The equilibrium energy density is evaluated 
from the thermodynamics of the dual Anti-de Sitter black hole, and is predicted to be 

e = ^-N^. (88) 

The shear viscosity is evaluated as the zero-frequency limit of the absorption cross-section 
of gravitons by the black three-brane (with gravitons polarized parallel to the brane) [12], 
and is predicted to be 

V = \NlT\ (89) 

Hence, in this limit the ratio ^ = rj/w is predicted to equal l/(4irT). The diffusion constant 
and susceptibility for SU(A) R charges are extracted from the pole structure of the thermal 
correlation function of SU(A) R charge densities [25], and are predicted to be 

A = ± (90) 

and 

Xn = \N 2 c T\ (91) 

The supercharge diffusion constant D s has, to our knowledge, not yet been evaluated, but 
should be computable by the methods of Ref. [25] . 

Although these strong coupling (and large N c ) limits of transport coefficients, suscepti- 
bility, and energy density cannot be compared with direct field-theoretic calculations, they 
may be used as input into the effective hydrodynamic theory (whose validity is independent 
of coupling and N c ). Inserting the AdS/CFT predictions (88-89) into the result (46) yields 
the completely explicit form 



l - J dx° d 3 x e^ x ° ([T l3 (x),T kl (0)]) = (l5 lk 5 3l + l5 u 5, k -l5 t ,5 kl )T 3 

1 (92) 



7r 2 56 + 6 3/2 tt\uj\ 

for the spectral density of the stress-stress correlator in the strong coupling, large N c limit. 
This correlator is a non-trivial probe of the real-time dynamics of strongly-coupled super- 
symmetric Yang-Mills plasma. Because the AdS/CFT prediction of the supercharge dif- 
fusion constant D s is currently lacking, we cannot give an equally explicit result for the 
i?-symmetry current spectral density (86). 
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The result (92) displays the expected 1/N% relative suppression of the non-analytic part, 
and also shows that the ratio of the leading analytic to the leading non-analytic piece does 
not depend on the coupling constant Qym of the gauge theory in the strong coupling limit. 
In the language of the AdS/CFT correspondence, a' /R\ dS = 1/ a/ Qym-Nc , and <iirg s = g Y M- 
The ten-dimensional gravitational constant k scales as g s a' 2 ~ R\ds/N c - Thus, it is 
string/supergravity loops, proportional to k 2 , which are responsible for 1/N 2 corrections 
in field theory correlators. The corresponding one-loop amplitudes are not straightforward 
to evaluate, given that even the scalar propagator on the plane-symmetric AdS black hole 
background is not known analytically. 37 However, it would be quite interesting to reproduce 
hydrodynamic results like Eq. (92) directly from supergravity loop corrections. Just veri- 
fying the form of the small-frequency non-analyticity from the gravity amplitudes, without 
computing the precise coefficient, would be a worthwhile goal. 

Finally, given that in various AdS/CFT-like scenarios, finite-temperature field theories 
are believed to be exactly equivalent to string theories on backgrounds with thermal horizons, 
it is natural to expect that there must exist a correspondence between appropriate effective 
theories as well. Namely, string theories on thermal backgrounds should have low-energy 
descriptions which are dual to effective long distance descriptions of finite-temperature field 
theories, specifically hydrodynamics and kinetic theory (in weakly coupled regimes). As 
we have seen, constructing a string theory dual of hydrodynamics will necessarily require 
incorporating string loop effects in order to reproduce aspects of long distance dynamics, 
such as long-time tails, which are sensitive to the non-linearities of hydrodynamics. 
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APPENDIX A: LONG-TIME TAIL IN STRESS-STRESS CORRELATOR 

In a theory, which is not necessarily scale- invariant, inserting the stress tensor constitutive 
relation (21) into the zero-wavenumber connected stress-stress correlator gives 

+ A & H £ B ££ ^ q (t) + m m 5 kl B„, ££ (t) 

+ j ^ S kl B ££ , ££ (t) + S* 5 kl B PaPbiPcPd (t) , (Al) 

37 For progress in this direction, see Ref. [26] . 
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where ffij = lXi Xj i I Xi Xj _ ? ,2 Xij X anf ] 
Wlieie n mn — 2 u m u n ' 2 U n u m u s u u mn, <AI1U 



Bir^i ,ee\t) 



d 3 X (±{TT%x)tf(t,x), 7T fc (0, 0)^(0, 0)}) c 

d 3 x (±{*%x)*S(t,x), £(0,0) 2 }> conn , 



5, 



,^{t) = J d 3 x {\{e{t,xf, 7r J (0, 0)^(0, 0)}> c 



B 



D :: . :: (f) = I <f'x n{v.xf. s(o,o) 2 }> conn . 

d 3 a; (| {p a (t, x) p b (t, x), p c (0, 0) p d (0, 0)}) £ 



Papb,PcPd 



it) 



(A2) 
(A3) 
(A4) 
(A5) 
(A6) 



Evaluating these correlators in the hydro dynamic regime, as described in section II E, yields 

d 3 k 



(27T) 



[A 7I . l7r fc(t, k) A^j^iit, —k) + A^i^i^t, k) A n j n k(t, —k 



w 2 T 2 [5 ij 5 kl + 5 ik 5 jl + 5 il 5 jk 2 <^'<5 fc ' + 7 (<W + W fe ) 
15 ~ 



(87T 7r; |t|)3/2 

/d :i k 
(27r) 3 ^ ^ 7ri£ ^' ^)^ 7rJ e(^' ~~ 



(47T 7s |t|) 3 / 2 
d 3 fc 



(A7) 



w 2 T 2 

3v 2 (47T 7s |t|)3/2 ' 
d 3 k 



B £ e,ee(t) 



(A8) 



(27T) 
2 T^2 



2A ee (t, fc) A ee (t,-fc) 



w 2 T 



1 



5 



PaPb,PcPd 



(t) 



vj (4vr 7s |t|) 3 / 2 ' 
d 3 k 



(A9) 



(27r) 3 ^ papc ^' ^ ^Pbpd(t, —k) + A PaPd (t, k) A PbPc (t, —k)] 

1/2 1/2 1/2 1/2 ba'd 8b>d' + ^a'cZ' <W 
Xaa' Xb&' X CC ' [ 47r ( jDa , +j D 6 ,)|t|]3/2 " 



(A10) 



In writing the result (A10), we have chosen to use a basis for conserved charges in which 
the symmetric matrix x~ 1 ^ 2 Dx 1 ^ 2 is diagonal and has eigenvalues {D a }. 
Putting everything together yields 

v-'(HT«(i),T"(o)}> = v«r^ + iM»E 



2T 2 



15 

T 2 



(47r 7s |t|)3/2 (87r7,|t|) 3 / 2 



(1 - 3v 2 + |u B - 2 ^) 2 4(l-3u, 



2\2 



(47T 7s |t|)3/2 

up to terms vanishing faster than |t|~ 3 / 2 as \t\ — > 00. 



+ 



(87^1)3/2 



(All) 
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